function derivs=Phi_d2(vars,params,index)
c_ba1=vars(1);
muz_cu=vars(2);
d_cu=vars(3);
n_cu=vars(4);
paistar_cu=vars(5);
a_cu=vars(6);
sigma_perturbation=vars(7);
BETTA=params(1);
B=params(2);
CHI=params(3);
CHI0=params(4);
THETA=params(5);
DELTA=params(6);
ALFA=params(7);
PHI=params(8);
PHIzero=params(9);
ZI=params(10);
KAPAw=params(11);
ETA=params(12);
RHOR=params(13);
PHIpai=params(14);
PHIpai_1=params(15);
PHIy=params(16);
PHIy_1=params(17);
PHIc=params(18);
PHIc_1=params(19);
PHIl=params(20);
NU=params(21);
U0=params(22);
U0d=params(23);
RHOz=params(24);
RHOd=params(25);
RHOn=params(26);
RHOp=params(27);
RHOa=params(28);
Kss=params(29);
OUTPUTss=params(30);
PAIss=params(31);
MUZss=params(32);
Css=params(33);
AA=params(34);
lss=params(35);
Rss=params(36);
Wss=params(37);
full_rows=zeros(0,1);
full_cols=zeros(0,1);
full_vals=zeros(0,1);
derivs=sparse(full_rows,full_cols,full_vals,0,49);
